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^ \ Abstract 

We investigate the gravitational energy-momentum distribution in 
the space-time of two black holes in circular orbit, in the context of 
the teleparallel equivalent of general relativity. This field configuration 
is important because gravitational waves are expected to be emitted 
in the final stages of inspiral and merger of binary black holes. We 
address an approximate solution of Einstein's field equations that de- 
scribes two non-spinning black holes that circle each other in the xy 
plane, obtain the total energy of the space-time and verify that the 
gravitational binding energy is negative. We show that gravitational 
radiation is emitted as long as the separation between the holes de- 
creases in time. If the black holes are spinning and circle each other, 
it has been found in the literature that, during the pre-merger inspi- 
ral, they bob up and down sinusoidally. The understanding of this 
phenomenon requires the understanding of the gravitational energy- 
momentum flow in the space-time of binary black holes. For the time 
dependent metric tensor of a general binary black hole system, the 
non-vanishing of the gravitational momentum may explain the bob- 
bing of spinning black holes. 
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1 Introduction 



The physics of black holes coalescence is presently being intensively investi- 
gated. So far it is not known any exact solution of Einstein's field equations 
that describes the inspiral and merger of black holes. The progress in this 
field is due to post-Newtonian approximation methods and to nonlinear nu- 
merical simulations of the evolution of binary black hole space-times [U [2], [3] . 
It is expected that binary black hole mergers will provide important informa- 
tion about the strong-field, nonlinear nature of the gravitational field, and 
will become a promising source of gravitational waves. Effects like radiation 
of mass (energy), linear and angular momentum are likely to take place in 
the final stages of black hole mergers. In particular, radiation of linear mo- 
mentum by a binary black hole is related to the recoil of the final remnant 
hole (see refs. [U EJ [6j [7J and references therein). If the remnant black hole 
acquires linear momentum, the latter should cancel with the linear momen- 
tum of the field in order to comply with the conservation of the total linear 
momentum of the configuration. 

An intriguing and interesting phenomenon is the orbital motion of two 
identical, spinning black holes in quasi-circular orbit, with oppositely directed 
spins restricted to the orbital plane [5J [U [9] . During the pre-merger inspiral 
the orbital plane of the binary black holes (plane xy, say) carry out a move- 
ment up and down along the z direction, i.e., the two black holes bob up and 
down sinusoidally and synchronously. After the merger the remnant black 
hole acquires a recoil velocity, and in realistic situations it may be ejected 
from the nucleous of the host galaxy jU [5j [7J . 

The conservation of the total linear momentum of the space-time implies 
that the linear momentum of the field must be the same in magnitude, but 
opposite in sign, to the linear momentum of the binary/merged black holes 
[8]. Therefore the description of the linear momentum of the gravitational 
field in the space-time of binary black holes is important for the understand- 
ing of the physics of this configuration. Ideally one would like to know the 
details of the momentum flow between the fields and holes. These issues are 
mathematically intricate, but in principle they can be addressed (at least 
formally) in the framework of the teleparallel equivalent of general relativity 
(TEGR), provided a realistic post-Newtonian expression for the metric ten- 
sor is known. In the TEGR the expressions for the energy-momentum and 
angular momentum of the gravitational field are invariant under transfor- 
mations of the coordinates of the three-dimensional spacelike surface, under 



2 



time reparametrizations, but depend on the frame of an observer. But nor- 
mally the observer is stationary in the asymptotically flat space-time. The 
energy-momentum of any physical system in special relativity depends on 
the frame of an observer, and there is no special reason for dropping this 
feature when considering general relativity. The gravitational energy is the 
zero component of the gravitational energy-momentum four- vector, and thus 
it has standard transformation properties. 

In this paper, we first address an approximate solution [10J [TTJ of Ein- 
stein's equations that describes two nonspinning black holes in circular mo- 
tion, and evaluate the gravitational energy-momentum of the space-time. In 
this approximate model the orbital motion of the holes is restricted to the xy 
plane and the separation b between the holes is considered fixed (in the con- 
text of ref. [10J). For two black holes with individual rest masses mi and rri2, 
we will find that the average value (in time) of the total energy cP^ of the 
combined system in orbital motion is less than c(P 1 ^°' ) + P^) = {mi + m2)c 2 . 
The binding energy is negative, its expression is 

very simple and is in agreement with previous analyses. Assuming that the 
separation b between the two black holes decreases in time, i.e., b < 0, the 
binding energy yields a positive flux of gravitational radiation. For a fixed 
value of b we arrive at a simple expression for the total flux of gravitational 
radiation. We will show that the average value of the radiation over a period 
is zero. The conclusion is that effective gravitational radiation takes place 
provided b 0. 

We also consider the general form of the metric tensor that describes the 
space-time of two spinning black holes in quasi-circular motion in the xy plane 
[HI [12], and evaluate the gravitational energy- momentum of the space-time. 
We obtain formal expressions for the linear momentum of the field along the 
x, y and z directions, contained within a large rectangular volume with sides 
a. This length is supposed to be much larger than the separation b between 
the two black holes. We find that, in general, the linear momenta vary with 
time. The dependence in time of the z component of the linear momentum 
is likely to be related to the bobbing of realistic binary black holes, prior 
to the merger. In our opinion, the present approach is better suited for 
this analysis, rather than the one based on pseudotensors [8] . Pseudotensors 
are quantities that depend on the choice of the coordinates of the three- 
dimensional space, and therefore they are not well behaved under coordinate 
transformations. The analysis developed in ref. [8] makes use of the Landau- 
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Lifshitz pseudotensor [13]. However, we have never seen a justification as to 
why one pseudotensor is better than another one. 

This paper is organized as follows. In section 2, we review the formula- 
tion of the TEGR, and show how the definition of the gravitational energy- 
momentum arises out of the field equations of the theory. We also present the 
expression for the fluxes of gravitational radiation and radiation of matter 
fields. In section 3, we describe the approximate solution for two black holes 
in circular orbit, and in section 4 we evaluate the energy and momentum 
of the gravitational field. We find simple expressions for the total gravita- 
tional energy, for the binding energy and for the total flux of gravitational 
radiation. In section 5, we consider the general form of the metric tensor for 
the inspiral of two spinning black holes that circle each other, and obtain 
the formal expressions for the gravitational momenta along the three spa- 
tial directions. Assuming the standard asymptotic behaviour of the metric 
tensor components, it will be clear that the momentum components of the 
gravitational field are time dependent, a fact that very likely explains the 
bobbing of the black holes. Finally we present our conclusions in section 6. 

Notation: space-time indices /i, u, ... and SO(3,l) indices a, b, ... run from 
to 3. Time and space indices are indicated according to \i — 0, i, a = (0), (i). 
The tetrad field is denoted e a M , and the torsion tensor reads T ativ = d^e av — 
d v e ati . The flat, Minkowski spacetime metric tensor raises and lowers tetrad 
indices and is fixed by r\ ah = e afl et >v g tlP = (— 1, +1, +1, +1). The determinant 
of the tetrad field is represented by e = det(e a M ). 



2 Energy- momentum in the TEGR 

In the teleparallel equivalent of general relativity the gravitational field is 
represented by the tetrad field e a M only, and the Lagrangian density is written 
in terms of the torsion tensor T aflv = d^e au — d u e a ^. This tensor is related 
to the antisymmetric part of the Weitzenbock connection T^ u = e aX d^e au . 
However, the dynamics of the gravitational field in the TEGR is essentially 
the same as in the usual metric formulation. The physics in both formulations 
is identically the same. 

Let us start with the torsion-free, Levi-Civita connection °uj^ ah) 
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^/jab 2^ ^iS^obc ^bac ^cab) j (1) 

^abc ^auifib^d^c ^c^d^Ch ) • 

The Christoffel symbols °r* and the Levi-Civita connection are identically 
related by 

= e a %e au + e aX (V^V 

In view of this expression an identity arises between the Levi-Civita connec- 
tion and the contorsion tensor K^ ab , 

°^fiab = , (2) 

where K^ ab = \e a x e b v (T XjXU + T vXjX + T^), and T XfXV = e a x T ailv . Making use 
of eq. (2) it follows that the scalar curvature R(e) may be identically written 

as 

eR(°u) = -e {^T abc T abc + l -T abc T bac - T a T a ) + 29„(eT") , (3) 

where e is the determinant of the tetrad field. Therefore in the framework of 
the TEGR the Lagrangian density for the gravitational and matter fields is 
defined by 



-ke(-T abc T abc + -T abc T bac - T a T a ) - -L 



1 , 

' M 



= -keZ abc T abc --L M , (4) 
where k = c 3 /16nG, T a = T b ba , T abc = e b ^e c u T a ^ v and 

j^abc _^rpabc _|_ rpbac rpcab^ _|_ ^^acrpb fj ^"] 10 '^ (5) 

Lm stands for the Lagrangian density for the matter fields. The Lagrangian 
density L is invariant under the global SO(3,l) group. The absence in the 
Lagrangian density of the divergence term on the right hand side of eq. (3) 
prevents the invariance of (4) under arbitrary local SO(3,l) transformations. 

The field equations derived from (4) are equivalent to Einstein's equa- 
tions. They read 
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e a \e blJ ,d u (eE 



b\u\ 



e(E bv a T hvil - \e ai [T hcd Y}> cd ) = ^eT a/ , , (6) 

where SLM/Se atl = eT afl . From now on we will make c — 1 — G. 

The definition of the gravitational energy-momentum may be established 
in the framework of the Lagrangian formulation defined by (4), according to 
the procedure of ref. [H]. Equation (6) may be rewritten as 



d u (eE 



a\u\ 



ee a M (t A " + T A ^) 



4k"" M 

where T x ^ = e a x T a>M and t x ^ is defined by 

t x » = k(4Z bcX T bc » - g x ^ bcd T bcd ) . 
In view of the antisymmetry property J] afJ-u = — J] ai/fl it follows that 



e e 



0. 



The equation above yields the continuity (or balance) equation, 



i- [ d 3 xee a Jt^ + T ^) 
at Jv 



dSj 



ee%(* 3 '" + T J >) 



(7) 



(8) 



(9) 



(10) 



Therefore we identify t x ^ as the gravitational energy-momentum tensor 

P a = f d 3 xee a Jt°" + T 0>M ) , (11) 
Jv 

as the total energy-momentum contained within a volume V of the three- 
dimensional space, 
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dSj (ee%F) 



;i2) 



as the gravitational energy-momentum flux [15], and 

^ m =ldSj (e e %n, (13) 

as the energy-momentum flux of matter. In view of (7) eq. (11) may be 
written as 



pa 



d 3 xd j U aj 



dSj U aj 



(14) 
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where IT ajf = — AkeH 0,0 ^. The expression above is the definition for the gravi- 
tational energy- momentum presented in ref. [IB], obtained in the framework 
of the vacuum field equations in Hamiltonian form. It is invariant under 
coordinate transformations of the three-dimensional space and under time 
reparametrizations. We note that (9) is a true energy-momentum conserva- 
tion equation. 

Finally we remark that in the absence of matter fields the total flux 
of gravitational radiation $ < -°- ) is related to the total gravitational energy 
according to = —p(°\ in view of eq. (10). 



3 The space-time of the binary black hole 

We consider the approximate metric for the binary black hole as described 
in ref. [10]. This solution was later re-analyzed and re-obtained in ref. [TT] 
by matching two perturbed Schwarzschild metrics to an asymptotically post- 
Newtonian construction for a binary black hole space-time. The two black 
holes have masses mi and m-i, and circle around each other in the plane xy. 
We restrict the analysis to the coordinates in the radiation zone, defined by 
r » Xqw [TO HB], where Xqw is the wavelength of the gravitational radia- 
tion. In the context of the present analysis the radiation zone is established 
in an equivalent way by r >> m 1; r » m 2 and r » b (see below the 
definition of b). We define 

m = mi + 7«2 , dm = mi — rri2 , U = . (15) 

m 

The circular, Newtonian trajectories of the black holes are 

n(f) = ^b(f), r 2 (t) = -— b(t) , (16) 
m m 

where 



b(t) = ri(t) - r 2 (t) = 6(coswt,sinwt,0) , (17) 

and 




is the orbital angular velocity. The separation b is defined by b = \r\ — r 2 
The velocity of the holes are 
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Vl 



~dt 



v 2 



dr 2 
~dt 



(18) 



from what we define 



v(t) = Vi - v 2 



smut, cosut, 0) . 



(19) 



In the radiation zone the metric components depend on b and v, which 
in turn depend not exactly on t, but on the retarded time r = t — r [in]. We 
further define 



m = m\ 1 — 



2b 



n 



(20) 



where r is a point of observation in space. 

With the help of the definitions above the metric tensor components for 
the binary black hole space-time in the radiation zone read pi)] 



goo 



9ij 



2m 2m 2 
■1 + - 



+ ^( 2 ( n ■ v) 2 - ^(n ■ b) 2 + ~(n ■ b)(n ■ v) + i[3(n • b) 2 - b 2 



+ 



u 5m [ , 

\ n 

r m I 



v) 



— (n.b) 2 -2(n.v) 2 - T 



|^(n.b) 2 -6(n.v) 2 -^ 



4/x 

9oi = 



+ -(n-b) 

r 

+^(n • v)[6 2 - 5(n • b) 2 ] + -i(n • h)[3b 2 - 5(n • b) 2 ] 



1 



(n • v) H — (n • b) 



S(n-b)^) 



b 3 



2u<5m/r. o 3m, 2 6. ... 
+— — 2(n- v) 2 - — (n-b) 2 + -(n-b)(n 

r m VI b 6 r 



+-[3(n-b) 2 -fe 2 ] \v 



+ 



4m, , . . . m 
■ — ( n -b)(n-v) + - 



b 2 



[n-bf 



. , 2m m 2 
Sij[l + — + 
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+ ~M n • v)2 - l? (n • b)2 + • b)(n • v) + ^ [3(n ' b)2 

+--{(- ■v)^(n-b) 2 -2(n.v) 2 + ^ 
r m I L o o . 




+-(n-b) -( n -b) 2 -2(n-v) 2 

+ -|(n • v)[6 2 - 5(n • b) 2 ] + l(n • b)[36 2 - 5(n • b) 2 ]}) 




(21) 



Our aim is to evaluate definitions (14) and (12) for the gravitational 
energy-momentum and the corresponding fluxes. These definitions are in- 
variant under global S0(3,l) transformations. Therefore they are frame de- 
pendent. However, P a is a vector under global Lorentz transformations. The 
frame dependence of the gravitational energy-momentum is understood by 
simply considering a black hole of mass m and an observer that is very distant 
from the black hole. The black hole will appear to this observer as a particle 
of mass m, with energy E = cP^ = mc 2 (m is the mass of the black hole 
in the frame where the black hole is at rest). If, however, the black hole is 
moving at velocity v with respect to the observer, then its total gravitational 
energy will be E = ^mc 2 , where 7 = (1 — v 2 j 'c 2 )^ 1 ^ 2 . This example is a 
consequence of the special theory of relativity, and shows clearly the frame 
dependence of the gravitational energy-momentum. The frame dependence is 
not restricted to observers at spacelike infinity. It holds for observers located 
everywhere in the three-dimensional space. 

In order to evaluate definitons (11-14) out of the metric tensor given 
by (21) we choose a configuration of tetrad fields that has a clear physical 
interpretation. In the framework of the TEGR the tetrad field describes both 
the gravitational field and the frame. For a given metric tensor there exists an 
infinity of possible frames, and each frame is characterized by six conditions 
on the tetrad field. Three conditions fix the kinematical state of the observer 
in the three-dimensional space (for instance, the observer may be stationary 
in space), and the other three conditions fix the orientation of the frame 
(alternatively, the frame may be characterized by the six components of the 
acceleration tensor 4>ab [1QJ ) . 
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Therefore tetrad fields are interpreted as reference frames adapted to 
preferred fields of observers in spacetime. This interpretation is possible by 
identifying the e(o) M components of the frame with the four-velocities of 
the observers, ^ = u^ 1 [19J. Here we will establish a set of tetrad fields 
adapted to static observers in spacetime. Thus we require em) 1 = 0. This 
condition fixes 3 components of the frame. The other three components are 
fixed by choosing an orientation of the frame in the three-dimensional space. 
Therefore e(o) ^ is parallel to the worldline of the observers, and e^k) ^ are 
the three unit vectors orthogonal to the timelike direction. We fix en-) M such 
that e(i) M , e( 2 ) M and e( 3 ) M in cartesian coordinates (and in the flat space-time 
limit) are unit vectors along the x, y and z directions. The tetrad field that 
satisfies these conditions is given by 



e a n(t,x,y,z) 



with the following definitions: 



[ABC 0\ 

D E F 

G H 

\ I J 



(22) 



.4 
B 

C 

D 

E 
F 

G 

H 

I 



[-g o) 1/2 , 

gat 
(-to) 1 / 2 

902 

(-goo) 1 / 2 

Aii 



-goo) 1 ' 2 
l 



x 2 

A 12 



-goo) 112 An 
-goo) 1/2 ^ 
l 



-flnnW 2 



An' 

yz _ A Vz 
A ll 



1/2 



'goo) 1 ' 2 g23^ 2 \ — gi?,^ 2 2 
An (AfiAi 2 -Af 2 )V2 



1 

A^ 



toAn - (-goo) (5i3 + 



(gggAu - £i 3 A 

(AllA2 2 



Af 2 ) 



nl/2 



(23) 
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The quantity \j is denned by A^- = goigoj — googij, and all metric compo- 
nents are obtained from (21). In the limit r — > oo the asymptotic quantities 
hoo, hn, h 2 2 and h 33 are defined by the expressions 



goo 


= -1 + /loo 


0n 




g22 


= l + h 22 , 


933 


= l + h 33 . 



(24) 



In terms of these quantities the asymptotic form of the tetrad field is reduced 
to 



e a n(t,x,y,z) 






V o 



-0oi 

x ~ 2 






-002 
012 



1 + 



fl22 









013 

023 
1 I hga 

~ 2 



(25) 



Expression (22) represents a frame that is adapted to static observers every- 
where in space-time. 



4 Gravitational energy of binary black holes 
in circular motion 

For a given space-time metric tensor and a given frame, the energy-momentum 
of the space-time is evaluated out of eq. (14). It reads 



P a = 4k 



(26) 



If the surface of integration S is fixed at spatial infinity, i.e., S — ?■ oo, P a yields 
the total energy-momentum of the space-time. The latter is the same for all 
tetrad fields that exhibit the same asymptotic behaviour. In particular, the 
energy-momentum obtained out of frames that are adapted to static observers 
at spacelike infinity coincides with the one obtained out of (22). 

Considering the tetrad field given by eq. (22), the gravitational energy of 
the space-time determined by eq. (21) is given by 
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= AkejdS ] e(e^\^ + e^ 1 T}^ + e^ 2 T, 2 ^) 

= 4k <f dS j e(ATP Qj + BE 10j + CS 2 ° 3 ) , (27) 

where A, £? and C are defined by (23), is calculated out of (5), and 

j = 1,2,3. We will evaluate the expression of P^ for a closed surface S in 
the asymptotic region r » m and r » b, which characterizes the radiation 
zone, and then we make S — > oo. The metric tensor may be decomposed 
as 9iiv — Vnv + hfj, v , and h^ v is of order 1/r at spacelike infinity. Thus we 
also have g^ v = rf v — where h^ u = rf^rf^h^. In view of the intricate 
structure of the metric tensor given by (21), in the evaluation of (27) we will 
make the approximation = r]^ a r] v ^r] Xl T la p 1 , making sure that eS" - 7 is of 
order (at least) 1/r 2 . After a large number of calculations, and taking into 
account the approximations explained above, we obtain the exact expression 
for the total energy of the space-time, 



p(0) 



(mi + m 2 ) 
r 8 fi 2 m 2 



1 - 



m x m 2 



3 b 2 



co cos(2wt) + 



2 (mi + m 2 )b. 
ix 1 '(Sm 



b 2 



muo\ 
~b~) 



cos(2wr) 



3 4 29 . 4 
+ - COS COT + — sm OUT 
8 8 



}si 



sm UT COS UT . 



(28) 



Assuming that the energy-momentum of matter fields vanishes for the binary 
black holes, then the expression above does represent the gravitational energy 
of the space-time. 

We will be interested in average values of time dependent quantities. 
Therefore we define 



<p(°) >=I FdrPW, 
T Jo 

where T = 2n /to. We easily obtain 

< p(°) >= ( mi + m 2 ) 



m x m 2 



(29) 



(30) 



2(mi + m 2 )b. 

We define the binding energy of the configuration according to =< p(°) > 
—(mi + m 2 ). We find 
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^ = -4^. (3D 

We note that Ej, is the the standard non-spinning part of the expression for 
the binding gravitational energy [201 El] for two black holes in circular orbit. 
This term is precisely the same as the first term in eq. (191) of ref. [22J. 

We evaluate now the total flux = — p(°) of gravitational radiation. 
After a large number of calculations we arrive at 



$(°) 



.2^2 



(bum 9 



3 b 2 w cos(2u;r) + ^ 



fi 2 (5m) 2 ( moo 2 



- cos(2ujt) H — cos 4 LOT 
2 8 



29 . 4 
+ — sin LOT 



| cos(2wt) 



2™2 



+ {--^(2u; 2 )sin(2u;r) + 



jj?{5m) 2 ( moo 2 



3 29 

- cos 3 lot sin lot H sin 3 lot cos o;r 

2 2 



b 2 \ b 
}si 



3 sin(2o;r) 



sinwT coswr . 



(32) 



An interesting consequence of the expression above is that the average value 
of $ < -°) over a complete cycle vanishes, 



<$(°)>=0. (33) 

We conclude that the orbital (stationary) motion of two black holes (or two 
point masses) on a plane produce the gravitational radiation given by eq. 
(32). However, the average value of this radiation vanishes. On the other 
hand, the situation changes if the separation distance b changes with time. 

In order to obtain (32) we have considered that b is not a function of t. 
This condition was assumed in ref. [ID] . According to ref. however, 
this condition may be relaxed. Thus we may admit that b varies slowly with 
time. Taking into account this assumption we obtain an expression for the 
total flux of the gravitational radiation. For arbitrary b(t) we have Co ^ 0. 
Consequently neither a complete cycle nor a period can be defined. In fact, 
the concept of average value does not apply to this case. $^°- ) is now clearly 
nonvanishing. We present here only the contribution to $^°- ) of the variation 
in time of the first term on the right hand side of eq. (28), namely, the 
variation in time of the binding energy, which does not even contribute to 
(32). We find 
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(34) 



This flux will be positive definite provided b < 0, which is the case for realistic 
binary black holes prior to the merger. 

We have evaluated the momentum component P( 3 \ which is the compo- 
nent oriented along the z direction. The integration is carried out over a 
finite cubic volume with sides a, such that a >> b. We found that P^ = 0. 
This result is expected since the black holes are restricted to the xy plane, 
and there is no flux of momentum along the z direction, in contrast to a 
general situation to be addressed in the next section. When the integration 
is carried out over the whole three-dimensional space, the total momentum 
of the space-time vanishes, i.e., P^' = for i = 1, 2, 3. 

5 Arbitrary time dependent metric tensor 

The standard form of the metric tensor for the space-time of the inspiral of 
two black holes is given by [23J [121 E] 

goo = -\ + 2V-2V 2 + 8X 
goi = -AVi-8Ri 

9ij = 5 lJ (l + 2V + 2V 2 )+AW lJ . (35) 

The form of the potentials V, Vi, X, and Wy may be obtained in the three 
references indicated above. Here we will just assume that V, Vi, Ri and X 
behave as 1/r at spacelike infinity, whereas W%j behaves as 1/r 2 . The expres- 
sion of the energy-momentum P a will be given in terms of these potentials. 
The use of the explicit form of the potentials yields a rather intricate form 
of the energy-momentum. The spinning nature of the solution is manifest in 
the potentials V, Vi and W„ [T2] . 

In cartesian coordinates a stationary observer in space-time is described 
by a frame very similar to eq. (22). In the metric tensor given by (21) we 
have #03 = 0, which is not the case for (35). The frame that (i) yields (35), 
(ii) is adapted to stationary observers in space-time, i.e., em) 1 = 0, and (iii) is 
oriented along the x, y, z directions at spacelike infinity, i.e., i(t, x, y, z) = 
b\ when r — > oo, is given by 
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e ,, = 



M 
o 



V o 



B C 

D E 

G 





J \ 

F 

H 
I ) 



(36) 



where the following relations are satisfied, 



A 2 = 


-goo 


AB = 


-001 


AC = 


~902 


AJ = 


~903 


-B 2 + D 2 = 


9n 


-BC + DE = 


912 


-BJ+FD = 


913 


^C 2 + E 2 + G 2 = 


922 


-CJ + EF +GH = 


923 


-J 2 + F 2 + H 2 + I 2 = 


933 ■ 



(37) 

The relations above allow to obtain all tetrad components in terms of the 
metric tensor components. 

We present below all components of the energy-momentum obtained out 
of the tetrad field (36), assuming the asymptotic behaviour of the potentials 
as explained above. We have discarded several terms that fall off as 0(l/r 3 ) 
or faster. For we obtain 



P(°) = -k{ lim / dydz [8,(922 + 9^)] 

/'OO 
dxdz [d 2 (gn + #33)] 
-00 

/OO 
dxdy [d 3 (g u + g 22 )] 
-00 

/OO 
dydz [Cd (g 12 ) + Jd (gn)\ 
-OO 

/OO 
dxdz [Bd (g 12 ) + Jd (g 23 )} 
-00 
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dxdy[Bd (g 13 ) + Cd (g 23 )} 



+ lim , 

Z^ZLOO J 

/oo roc 
dydzBdo(gn) + lim / dxdzCdo(g22) 
-oo y->±°oj_ 00 

oo 



+ Z ^ 00 J_ d xdyJd (g 33 )} , 



(38) 



where 



/oo f'OQ roo roo roo 

dydzF(x,y, z) = j dy j dzF(oo,y,z)— dy / dzF(—oo,y,z). 
-oo J—oo J—oo J—oo J—oo 



Equation (38) reproduces (28) if we reduce the metric tensor (35) to the 
form given by (21), in which case we make J = 0. 

In view of the asymptotic behaviour of the metric tensor (35), the total 
momentum of the space-time vanishes. In the expressions below for pW ) we 
formally integrate over a finite surface So of a large rectangular volume with 
sides (2xo, 2y , 2z ), such that x = y = z = a, and a is much larger than 
the separation of the black holes. We arrive at 



^ dy dziVidoV!) + / dzdx(V 1 d V 2 ) 

So J So 



+ 

'So 

+ 4k 



P (1) = -32k 

f dxdyiV.doVs 

JSo 

f dydzd (V + W 2 2 + W 33 ) 

J So 

/ dxdy(d Wi 3 )- dzdx(d Wi 2 ) 

J So J So 



(39) 



P (2) = -32A; 



So 



dydz(V 2 d V 1 )+ j dzdx(V 2 d V2) 

J So 



+ [ dxdy(V 2 d V 3 ) 

J So 



+ 4k 



So 



dzdxd (V + W 11 + W 33 ) 



dy dz(d W 12 ) - dxdy(d W 23 ) 

So J So 



(40) 
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p(3) 



where 



-32A; 



dydz{V 3 d V 1 )+ [ dzdx(V 3 d V 2 ) 
+ / dxdy(V 3 d V 3 

J So 

I dydzd {V + W n + W 22 ) 

J So 



I S 

+ Ak 



dzdx(d W 23 ) — / dydz(d Wi 3 ) 

So J So 



(41) 



•So 



dydzF(x, y, z) 



yo rzo ryo rzn 

dy dzF(x ,y,z)- dy dzF(-x ,y,z) 
-yo J—zq J—yo J—zo 



etc. An important conclusion that we can draw from the expressions above 
is that the potentials X and Ri do not contribute to the momenta. 

Considering the post- Newtonian potentials V, Vi and Wij, we find that the 
component of the gravitational momentum is, in general, nonvanishing 
for a finite volume of the three-dimensional space, and exhibits a dependence 
in time. Therefore, it yields a momentum flux (p^ = — p( 3 ) which, in turn, 
is likely to be related to the bobbing of the black holes. Unfortunately the 
post-Newtonian potentials of ref. [23J are not suitable for the present anal- 
ysis, because they are valid only in the near zone, and we are interested 
in the expressions of the potentials in the radiation zone. Some of the po- 
tentials presented in the latter reference diverge with the increasing of the 
radial distance r, a feature that prevents us from calculating all momentum 
components. 

The post-Newtonian potentials V, Vi and depend on the spins S x , S y 
and S z of the black holes, which are time dependent functions [12]. The 
explicit form of these functions, for given initial conditions, and an analytic, 
widely accepted form of the post-Newtonian potentials, are not available in 
the literature. For this reason, we cannot proceed and obtain the detailed 
form of expressions (39), (40) and (41). 
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6 Concluding remarks 



In this paper we have analyzed the metric tensor for the nonspinning black 
hole binary in circular orbit in the xy plane, in the context of the teleparallel 
equivalent of general relativity. This metric tensor is an approximate solution 
of Einstein's equations in which the distance between the holes is constant 
in time. 

We have also addressed the general post-Newtonian form of the metric 
tensor that describes the inspiral and merger of two spinning black holes. 
The total energy-momentum of the space-time may be expressed in a simple 
form in terms of the metric tensor components, and may be easily computed 
provided the post-Newtonian potentials are given. We have found that only 
the potentials V, Vi and Wij contribute to the momenta. It is very likely 
that the time dependence of the momentum component P^ (the momentum 
component oriented along the z direction) is related to the bobbing of the 
spinning black holes with oppositely directed spins restricted to the orbital 
plane. 

The calculations of the total gravitational energy, and the corresponding 
flux for the nonspinning black hole binary in circular orbit, yield a quite 
interesting result. The total energy P^ and the gravitational flux cj)^ are 
given by eqs. (28) and (32). The former yields the known result for the 
mean value of the binding energy of the configuration, whereas the average 
value of cj)^ in time vanishes, < >= 0. It means that for two black holes 
in circular orbit, the average value of the gravitational radiation is zero. 
This result is consistent with the stationary character of the space-time. 
A nontrivial emission of gravitational radiation must necessarily be related 
to a loss of energy-momentum of the source (as the loss of mass described 
by Bondi's radiating metric [21]). If, however, the separation between the 
holes decreases in time, as in an actual evolution of the black hole binary, a 
nonvanishing (definite positive) flux of gravitational radiation is emitted. 

This result is conceptually different from the conclusion drawn from Ed- 
dington's spinning rod |25J , which was reconsidered by other authors . In 
the framework of pseudotensor definitions and of linearized general relativity, 
the quadrupole formula was obtained. The latter relates the energy loss of 
the system with the variation in time of the mass-quadrupole of the source. 
This formula was used by Eddington to deduce the energy flux generated by 
a rod that spins in the xy plane with angular frequency u. Let / represent 
the moment of inertia of the rod, and G the gravitational constant. The total 
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energy flux is given by [221 H3] 



dE 
~dt 



32G7 2 w 6 
5c 5 



(42) 



In contrast to the approach that allows to deduce the formula above, we 
note that the procedure that led to eqs. (32) and (33) is based neither 
on pseudotensors nor on the linearized form of the theory. Equation (10) 
is a true tensorial quantity. It is valid for any coordinate system of the 
three-dimensional spacelike hypersurface, and for finite volume V and cor- 
responding surface S. We believe that the reconsideration of Eddington's 
spinning rod in the present context would also lead to an equation of the 
type < (j)^ >= 0, provided the angular frequency is constant. 



References 

[1] F. Pretorius, Phys. Rev. Lett. 95, 121101 (2005). 

[2] M. Campanelli, C. Lousto, P. Marronetti and Y. Zlochower, Phys. Rev. 
Lett. 96, 111101 (2006). 

[3] F. Pretorius, "Binary Black Hole Coalescence". larXiv:0710. 13381 

[4] M. Campanelli, C. Lousto, Y. Zlochower and D. Merritt, Astrophys. J. 
Lett. 659, L5 (2007). 

[5] M. Campanelli, C. Lousto, Y. Zlochower and D. Merritt, Phys. Rev. 
Lett. 98, 231102 (2007). 

[6] J. A. Gonzalez, U. Sperhake, B. Briigmann, M. Hannan and S. Husa, 
Phys. Rev. Lett. 98, 091101 (2007). 

[7] L. Blanchet, M. S. S. Qusailah and C. M. Will, Astrophys. J. 635, 508 
(2005). 

[8] D. Keppel, D. A. Nichols, Y. Chen and K. S. Thorne, Phys. Rev. D 80, 
124015 2009. 

[9] G. Lovelace, Y. Chen, M. Cohen, J. D. Kaplan, D. Keppel, K. D. 
Matthews, D. A. Nichols, M. A. Scheel and U. Sperhake, Phys. Rev. 
D 82, 064031 (2010). 



19 



[10] K. Alvi, Phys. Rev. D 61, 124013 (2000). 

[11] N. K. Johnson-McDaniel, N. Yunes, W. Tichy and B. J. Owen, Phys. 
Rev. D 80, 124039 (2009). 

[12] G. Faye, L. Blanchet and A. Buonanno, Phys. Rev. D 74, 104033, (2006). 

[13] L. D. Landau and E. M Lifshitz, "The Classical Theory of Fields" (Perg- 
amon Press, Oxford, 1980). 



[14] J. W. Maluf, Ann. Phys. (Berlin) 14 (2005) 723 |gr-qc/0504077] . 

[15] J. W. Maluf, F. F. Faria and K. H. Castello-Branco, Class. Quantum 
Grav. 20 (2003) 4683. 

[16] J. W. Maluf, J. F. da Rocha-Neto, T. M. L. Toribio and K. H. Castello- 
Branco, Phys. Rev. D 65 (2002) 124001. 

[17] N. Yunes. W. Tichy, B. J. Owen and B. Briigmann, Phys. Rev. D 74, 
104011 (2006). 

[18] N. Yunes and W. Tichy, Phys. Rev. D 74, 064013 (2006). 

[19] J. W. Maluf, F. F. Faria and S. C. Ulhoa, Class. Quantum Grav. 24, 
2743 (2007). 

[20] R. Wald, Phys. Rev. D 6, 406 (1972). 

[21] S. Dain, Phys. Rev. D 66, 084019 (2002). 

[22] L. Blanchet, Living Rev. Relativity 9, (2006), 4. 

[23] H. Tagoshi, A. Ohashi and B. J. Owen, Phys. Rev. D 63, 044006 (2001). 

[24] H. Bondi, M. G. J. van der Burg, A. W. K. Metzner, Proc. R. Soc. 
London Ser. A 269, 21 (1962). 

[25] A. S. Eddington, Proc. R. Soc. A 102, 268 (1922); "The Mathematical 
Theory of Relativity" pp. 248-252 (Cambridge Univ. Press, Cambridge, 
1965), originally published in 1923. 

[26] F. I. Cooperstock, Ann. Phys. (NY) 282, 115 (2000); F. I. Cooperstock 
and S. Tieu, Found. Phys. 33, 1033 (2003). 



20 



